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Abstract 

We investigate the effects of Yukawa couplings on the phenomenological predictions for a class of 
supersymmetric models which allows for the presence of complete SU (5) multiplets in addition to the 
Minimal Supersymmetric Standard Model spectrum. We develop a two-loop analytical approach to 
quantify the predictions for gauge unification including Yukawa couplings. The effects of the heavy 
thresholds of the model are also included. In some cases accurate predictions can be made for the 
unification scale, irrespective of the initial (unknown) Yukawa couplings, so long as perturbation 
theory remains valid. We also consider the limit of a large number of extra states and compute the 
predictions in a resummed perturbation series approach to show that the results are stable in this 
limit. Finally we consider the possibility of making predictions for the case the gauge and Yukawa 
couplings enter the non-perturbative domain below the unification scale and estimate the errors 
which affect these predictions. 
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1 Introduction 



In a recent paper |T|] the phenomenological implications of a class of string motivated supersymmetric 
models based on level- 1 heterotic string models with Wilson line breaking of the underlying Eg symmetry 
were investigated. These models (henceforth referred to as Extended Minimal Supersymmetric Models 
(EMSSM)) contain additional vector-like representations filling out complete five and ten dimensional 
representations of SU (5) even though the gauge group is just that of the Standard Model[|. This class of 
string models has several advantages for model building |^]; the Wilson line breaking necessary to break 
the E$ symmetry in level- 1 theories offers an elegant way out of the doublet-triplet splitting problem 
encountered in GUTs. It can also lead to the standard unification of the gauge couplings without the 
need of a Grand Unified Group below the compactification scaleQ. In this case the string prediction 
for the scale of unification of the gauge couplings may be directly compared with the value obtained 
continuing the gauge couplings up in energy providing the possibility of a quantitative test of string 
unification including gravity. 

The effect of the additional vector-like states on 03 (M z ) and the value of the unification scale was 
explored in detail in|l]]. It was found that, working to two- loop order in the gauge sector, the unification 
scale is systematically increased by a small factor (about 3 or less, a function of the value of the unified 
coupling and of the number n = (77,5 + 3nio)/2 of pairs f] of additional complete SU(5) representations), 
while the strong coupling was systematically increased, taking it further from the experimental value 
0], 03 (M z ) = 0.118 ± 0.003. The increase in the unification scale does not, in fact, take it closer to 
the weakly coupled heterotic string prediction ||] because the increase is correlated with an increase in 
the unified coupling which always happens when additional matter is added and the heterotic string 
prediction also increases with the unified coupling. As a result there is still a discrepancy^ of a factor of 
10-20. However, in this analysis the effects of Yukawa couplings on the running of the gauge couplings 
were not included and we seek to include them in this paper. An immediate difficulty arises because in 
theories with additional massive vectorlike states there may be many new Yukawa couplings involving 
these heavy states. We shall consider two basic possibilities which give an indication of the general 
possibilities and uncertainties. In the first we include only the Yukawa couplings present in the MSSM 
which are responsible for the third generation masses. In the second we include couplings between 
massive Higgs doublets and the light quarks and leptons. This model goes some way towards explaining 
the light quark mass matrix via mixing in the Higgs sector Q as well as explaining the masses of the 
third generation, in a "fixed-point-scenario" We also comment on how our results would be affected 
if further couplings involving massive states were included. 

The paper is organized as follows. In Section 2 we review the effect of heavy thresholds on the 
running coupling constants and we analyze the renormalisation group evolution (RGE) predictions for 
the gauge couplings with Yukawa couplings effects included. In section 3 we derive an analytic form 
for 03 (M z ) and the unification scale as well as the decoupling scale of the extra matter and make some 
numerical estimates. In section 4 we consider the case of large n. Finally, in Section 5 we consider 
the case of unification at strong coupling. We discuss the differences between the (quasi)-fixed point 
approach || and the two loop perturbative approach results. Our conclusions are presented in the last 
section. 



1 To be precise these models contain additional vector-like states I + I , where /, / denote complete representations of 
517(5). The components ip of the additional complete SU(5) representations transform under SU(3) and SU(2) groups 
as follows: for t[> — d c : (3, 1) , t[> — I : (1,2) for the case I is the 5 dimensional representation of SU(5) and ip = e c : 
(1,1) , V> = u c : (3, 1) , ip = q : (3,2) for the case / is the 10 dimensional representation of 517(5) 

2 Alternative constructions were developed which can lead to non-standard U(l) normalisation even in level-1 theory - 
see discussion in ref. js| 



s n 5 = N 5 + Ng and n 10 = N w + N-^ 
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Ways to reconcile this scale discrepancy have been reviewed by Dienes 
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2 Heavy thresholds' contribution from the integrated NSVZ beta 
function 



In this section we make use of the Novikov-Shifman-Vainshtein-Zakharov (NSVZ) beta function |l0| 
for the gauge couplings, which was derived using holomorphicity arguments and the instanton calculus. 
We integrate it in the presence of the additional heavy states of the EMSSM model and then use its 
two-loop approximation to perform further calculations. 

The integrated expression for the gauge couplings' running, in the absence of these heavy thresholds 
due to the presence of the extra heavy states we consider, was presented in ref. 11] and was deduced 
without actually integrating the beta function, but by using physical arguments and the properties of 
the holomorphic gauge coupling; from that, the integrated form in the presence of our additional heavy 
thresholds can be easily "guessed" . For clarity, we prefer to re-derive it here by direct integration. For 

ive§ 

1-T(G>£)" (D 



a fuller discussion the reader is referred to the literature ]10| , 11, 12]. We have 
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d(ln /i) 2ir 
with the definition (p, is the running scale) 

NSVZ _ d\nZ a 
dliifi 

and where T(G) and T(R a ) represent the Dynkin index for the adjoint representation and for R a rep- 
resentation respectively (not necessarily the fundamental one), whose values are given in the Appendix. 
The above sum runs over all matter fields a in representation R a and this includes, for example, the 
extra heavy states (which we called ip -see footnote (1)), in addition to the low energy spectrum of 
MSSM. This gives 

fcr^Ll + T(G)§^ = 3T(G) - £ T(R a )(l - (2) 
d In fjL dm n 

which can be integrated from the high scale M down to the low scale fi with fi > Asusy to give 



2K{a-\M) - a' 1 + T(G) In ^ = 3T(G) In — - £ F T{R a ){\ - ^ Svz )d(lnfi) (3) 

a(n) n ah 

In the above equation a stands for any gauge group coupling. Any extra heavy state with mass ^(/x^) 
larger than /j,, decouples at its physical mass scale in the running of the canonical gauge coupling |l0| , [l2| 
and therefore, below scale, 7^ = 0. From eq.@ we obtain the result 



1/ n -3T(G) , M 

' 2vr / a(M) \ / 



1 ^ M 1 ^ Z^M) 

+ ^EW^4EW^ (4) 
where <j>'s stand for a MSSM-like spectrum. At this point we make the observation that 
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which is the mass renormahsation equation, where fjfL represents the bare mass of the SU(5) component 
tp. For the case of a Grand Unified Group Model, the gauge invariance principle requires the mass terms 
fJ^iptp be invariant, which is possible only if all bare masses are equal to a common value, \i g . We also 
consider, without any loss of generality that Z^(M g ) = Z^Mg) = 1 where M g is the value of unification 
scale. Taking M = M g , we find 



n 



1 u = oTVMJH —In ^— r7T+> — ^-^-ln 7— H In — 9 - (5) 

with n = (ns + 3rtio)/2, re 5 = ^5 + N§ and nio = iVio + -^Vjn- The last term in eq.(|5|) could stand for the 
Higgs in the adjoint representation or for the SU(3) Higgs triplet^. This formula is valid to all orders 
in perturbation theory, as long as Supersymmetry is not broken. This equation was used, in this form, 
as a check of the intermediate results obtained in [jjj . 

In two-loop order for gauge couplings running, Z factors depend on the mean mass of the extra 
heavy states which can be considered to be their common bare mass [i g since the difference would be 
of three loop order. Hence, in two-loop order /i g is the only mass scale associated with the heavy states 
in eq.(|5|). Therefore one sees that the decoupling of the heavy states happens, in two- loop order, at n g 
and not at each of the physical masses of the component fields of the massive SU(5) representations. 

3 Analytical Results from Renormalisation Group Evolution 

In this section we consider the effects of the Yukawa couplings on the running of the gauge couplings 
and their effects on the predictions for the unification scale, the strong coupling at electroweak scale 
and the scale where the heavy states decouple. In the case the Yukawa couplings are negligible, we 
recover previous two- loop analytical results M. 

To keep our approach general and to exhibit the effects Yukawa couplings have on the running of 
the gauge couplings, it is convenient to use eq. (|5j) for the running of the gauge couplings. We should 
mention that in two-loop order this formula has no regularization ambiguities, which only arise at three 
loop and beyond [13|. The advantage of eq.(||) for the evolution of the couplings is that it expresses the 



running of the gauge couplings in terms of the wavefunction renormalisation coefficients Z^ . Unlike the 
gauge wave function contribution which is one-loop exact, one still needs to use perturbation theory to 
compute the Z^ order by order for the matter wave function. However, for a two loop approximation 
for the gauge couplings, only a one-loop calculation of Z^ is required, simplifying the calculation. Their 
evolution from the unification scale where they are normalised to unityQ, to the decoupling scale \x g of 
the heavy states is affected by the presence of n = (715 + 3nio)/2 pairs of complete SU(b) multiplets 
which affects the anomalous dimensions of the matter fields]]. 

Below the supersymmetry breaking scale we have to include the effect of the low energy supersym- 
metric thresholds. Their effect at one-loop, which corresponds to a two-loop effect if the splitting of the 
super-partner masses is of radiative origin, we denote by <5j. This term is needed because eq. (|5j) is valid 
only as long as supersymmetry is unbroken []. 

5 These exotic Higgs effects will not be considered further in this paper 

depend on the unification scale! 
7 The running of the gauge couplings depends on ns + 3nio only and not separately on ns or nio (l[ . 
s Si include regularisation scheme conversion factors as well. 
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With these considerations we rewrite eq. (JsJ) for the three gauge couplings (j is a generation index) 
in the following way: 
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where 6i = 33/5, b 2 = 1, 63 = —3 and where /9j ^. = T(R\), i = {1, 2, 3}, are the contributions to one- 

loop beta function^of the matter fields <pj (j=generation index), while Pi t9 = T l (G) is the one-loop beta 
function for the pure gauge (+gaugino) sector; the Higgs (+higgsino) sector contribution is included 
separately via the terms proportional to A,Hi,2 considered. Thus we have 



/ _3_ j_ 3 4 1 \ 
'10 10 5 5 5 \ 



\ f 



V 1 



1 1 



Ji 



i,<t>j 



( \ 

-6 
V "9 J 



/ in \ 



(7) 



independent of the values of j. The field <f>j runs over the set 
with j as generation index. 



V J 

{h,qL,eR,u R ,d R }j, in this order, 



3.1 a 3 (M z ), M g and fi g . 

To compute the two-loop values for the strong coupling at electroweak scale, the unification scale and 
the decoupling scale n g of the extra massive states^, we use the two step method developed in reference 

i- 

The first step eliminates the low energy super symmetric thresholds dependence by expressing our 
results as a change to MSSM predictions for the strong coupling and the unification scale; this is possible 
because, to two loop order for the RGE, the <Vs are the same in both cases. Note that this means that 
our results will be expressed in terms of a°, M° and a°{M z ) which all have the 5i dependence included. 

The second step is more technical and consists of replacing the arguments of the two-loop log's by 
the one-loop approximation. This will, of course, generate the evolution of the couplings correct to 
two-loop order. 

To implement the first step note that in the MSSM we have 

<-\M z ) = -5 l + a ^ + ^ln^-^lnZ^(M z )-^lnZ° H2 (M z ) 

i =1 <f>, 

where we employed the index "o" to label all MSSM related quantities, to distinguish them from those 
of the extended model which introduces additional heavy states. Note that Z° are normalised to unity 
at the M;, Z°{M° g ) = Z°(0) = 1. 

9 We also used that the one-loop beta function is b = — 3T(G) + z \2T(R ( p), where the sum runs over all chiral super- 
multiplets in representation R^. 

10 fig is the decoupling scale only in two loop order 
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As is commonly done in all "bottom-up" approaches, we impose the conditions a±(M z ) = a\{M z ) 
and a2(M z ) = a 2 (M z ) and equal to their experimental values. We now determine the change of the 
strong coupling at electroweak scale; from @, (|8|) we get that (i=l,2): 
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(10) 

-l/a§(M 2 ). 



where by 5a^ 1 (M z ) we denoted the (two-loop induced) difference Sa^ (M z 
We note that the term ln(Qg(M 2 )/a3(M 2 ) can simply be neglected in the two-loop approximation as 
ajj(M 2 ) and a 3 (M z ) are equal in one-loop, and therefore this term would bring a higher order correction. 
We solve the system of eqs. (^), ( JlOj ) for 5a^ 1 (M z ), M g and n g in function of a g to get: 
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where <r^. = Pi^fa - b 3 )/(bi - b 2 ) + ^2,^(&3 - bi)/(h - b 2 ) + /33,<^ and we have a similar definition 
for crflj, crfl- 2 and <r 5 . Thus 
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(12) 



We observe that is negative for SU(2) doublets and positive for SU(2) singlets and these will 
therefore have opposite effects on the values of a 3 (M z ), for same sign of the associated logarithmic 
factor. Moreover, the contribution to the strong coupling depends on whether Z's are larger or smaller 
than unity, which is determined by the relative effect of gauge and Yukawa contributions to their 
running. This should become clearer later. The overall effect on the strong coupling depends however 
on the relative magnitude of the wavefunction coefficients, Z, of the extended model compared to those 
of MSSM, Z°. 

Using the notation A/3i 2| ^ = /3i,<^ - 02,^ and A/3i2 t H 1>2 = Pi,H 1>2 ~ fa,Hi,2 tne change in the 
unification scale factor is given by (see eqs.(S), fllO|)): 
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Note that A/3 12j( / )i has negative sign for SU(2) doublets and positive for 577(2) singlets (just like ) 
and hence, for same sign of the logarithmic factor, these terms will drive the unification scale in opposite 
directions. Moreover, comparing eqs.(|TTD and ( |13|), we observe that the unification scale M g is increased 
for positive A/?i2 eq.(13), while alpha strong eq.(|ll|) is also increased^ for a given (positive) sign of the 



1 A/3i2,0 has same sign as cr^. 
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log's. Thus, we already see it is difficult to decrease the strong coupling and increase unification scale 
simultaneously. One could argue that Z factors have an implicit unification scale dependence as well, 
and that our above explanation does not apply; however this dependence comes in under the log, and 
is therefore very mild. Note that this applies independent of the type of Yukawa interaction, implicitly 
present in the values of Z coefficients. This seems to be generic to models which consider complete 
additional representations and it was also explored in Q in the absence of Yukawa effects with same 
conclusions. The above discussion should become clearer later when a more quantitative study will be 
made. 

Finally, the decoupling scale [i g is given by (see eqs.(^), (10)): 
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where 
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with F = {(pj, Hi, H 2 } and b[ = bi + n. 

To simplify our expressions further we need to compute the wavefunction renormalisation coefficients. 
Above the mass of the additional states in the extended model the running of these coefficients is 
altered implying that the anomalous dimensions of the matter fields are changed. We have (with 
t = l/(2vr)ln(sca/e)) 



-±TT lF = 2j2C k (F)a k (t)+ £ A v {F)y v {t) + 0{o? 



(16) 



k=l 



u=T,b,t 



In the above equation Cj(F) is the second order Casimir operator for the field F for the group SU(j) 
and is generation independent. The coefficients A U (F) depend on the type of the superpotential of the 
model. For various cases we consider, their values are given in the Appendix. 

The above equation can easily be integrated to give full analytical expressions for the wavefunction 
coefficients, valid in one- loop approximation which is consistent with a two-loop running for the gauge 
couplings. The solution is of the type: 



Z F (t) = Z$(t) x Z F (t) 



(17) 



where Zp and Z^f are given by the gauge and Yukawa running parts respectively, with Zp(0) = 1 and 
Zp(0) = 1. Zp is determined by the gauge group; for the extended model we get from eq.([l6|) and ([l7]) 
that 



Z$(M Z 




2g fc (F) 



(18) 



Note that for one-loop running of the coefficients Z (consistent with two loop approximation for gauge 
couplings), the decoupling scale of the heavy states ip can be considered to be \i g rather than a physical 
mean mass of the multiplet, as the difference is a radiative effect and therefore represents a higher 
(two-loop) correction to Z's or three loop order for the gauge couplings. This justifies the last equation 
we wrote above. 
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Similarly, we have that in MSSM, 

2C fc (F) 

(19) 

The coefficients Zp can also be expressed in a form similar to that of Zp (eq. QT8|)), in terms of 
Yukawa and gauge couplings at initial and final scale only; to do this one only needs the one-loop 
running for Yukawa couplings; from this, any Yukawa coupling can be expressed (see Appendix) in 
terms of the derivative of \xiy T fi,t and derivative of lncti where the latter is obtained by replacing ati 
with (1 / bi) d In ati/dt, consistent with one-loop running for Yukawa; hence, the integral / y u dt can be 
performed analytically and therefore the differential equation of Zp's (similar to (|i~6| ) without the gauge 
term) can be integrated to give ZpS in one-loop approximation. 

A final comment is in order here: we see that the powers entering eq.(|l8l) depend on n as n/{bi + n). 
This means that in the limit of large n the contribution of these terms to the running of the couplings 
brings in only a small depedence on n. We will later see that a g /a,k([i g ) itself is also relatively stable 
with respect to n. Next, from the running of the gauge couplings eqs. (||) and @ we also observe 
that the explicit n dependence comes in under same form for all the couplings as n/(27r) ln(M g /fj, g ), 
which means that the relative behaviour of the gauge couplings evolution is affected only through Z 
factors. With a\{M z ) and 02(^2) fixed to their experimental values, and the above observations we 
conclude that the predictions for the strong coupling at electroweak scale will be stable for large n. 
This observation applies to other predictions we will make as well, such as the unification scale and 
the bare mass, [i g . For the absence of Yukawa couplings' effects this behaviour was shown in JjJ]. We 
will see that this remains true in their presence too, although a dependence of the particular type of 
superpotential might be expected. 
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3.2 Analysis for 3rd generation Yukawa couplings only. 

We will consider first the case when the superpotential is the same in both MSSM and in the extended 
model. In eqs. (|TT|) , (|l3[), (|i~4|), we replace the gauge part of Z coefficients, Zp, with their one-loop 
analytical expressions, consistent with two-loop running for the gauge couplings. For the Yukawa part 
of Z coefficients, Zp, from their evolution equation 

j t \nZ F '(t)= £ A„{F)y v (t) + 0{a 2 ) (20) 

u=r,b,t 

after using the relations among the coefficients A U (F), (not affected by the presence of the extra 
matter 12 ), one can get the following relations, valid at any scale between M g and M z 
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where we used, without any loss of generality, the convention Z (0) = 1 for any field. Making use of 
these relations, we get from eq.(|iT|) that 

5a^(M z ) = -^ln 

12 Their expressions are presented in Appendix 
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where Z Y , Z oY factors are evaluated at M z scale and where we defined 
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ln(a 9 /<x,(/x 9 )) by 



-1, 



(M 2 ) — ag 1 (M Z ). To get the above formula we replaced the terms of the form 



In 



aj(ng) 



In 



M n 



l + a -Mn — y - 



'2tt 



In 



b' : 



1 + ^ 



a" 



(24) 



where we made use of the fact that a g 1 
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+ n/ (27r) ln(M g /[ig) = in one-loop order |^| which is 
obtained from equations (^) and (10) after noticing that, in one-loop approximation, all wavefunction 
renormalisation coefficients are equal to 1. 

Similar results hold for the unification scale and for the decoupling scale of the extra matter. We 
have from eqs.(|l~3l), (|2~T|) 
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where Z Y , Z oY factors are evaluated at M z scale and 
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Finally from eq.fllil) and (^| 
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In the above expressions eqs.(^), (^), (p7|), Z y factors are evaluated at the electroweak scale, and 
normalised to unity at the unification scale M g and the same is true for the Z° coefficients^ If we 
set all "Yukawa" wavefunction coefficients Z^f and Z F Y equal to unity, we recover the results Q we 
obtained previously in the absence of any Yukawa couplings^. As the coefficients Zp > 1 (unlike 

13 In fact a" 1 -a^ + &i/(2ir) In [M g /M°) + n/(2n) In (M g /(j, g ) = in one loop and we further have In [M g /M°) = 
in one-loop because the change of unification scale is a two-loop effect. 

14 We also have that Sa^ 1 is non-zero in two-loop only, being in one-loop. 
15 These are normalised to unity at Mg ! 

16 In ref. jl[ the results of the extended model were compared with those of MSSM without top, bottom, tau Yukawa 
effects. 
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Zp < 1) (see eqs.(|l6|), (pO|)), we note that their effect is to decrease the strong coupling at electroweak 
scale. However, we must evaluate the relative effect of these coefficients Z Y to those of MSSM, Z oY \ this 
will be done below, by expressing them in terms of Yukawa couplings, assuming the same low energy 
input for Yukawa couplings in both EMSSM and MSSM models. This means we take same values at 
the electroweak scale for top, bottom and r Yukawa couplings in both models. Using the expressions 
of Z Y and Z oY , the result will be expressed as a function of ratios of the Yukawa couplings, following 
the details presented at the end of the previous subsection, after eq. (|i~9|) . 

For the strong coupling eq.(22) we obtain the following result: l/a^(M z ) = l/a%{M z ) + 8a% {M z ) with 
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j = {1,2,3}. The terms which involve a g /a° in the r.h.s. of eq,(|30"D drive a^{M z ) above MSSM value. 
We also see that, imposing the same electroweak scale values for Yukawa couplings in both models, the 
overall effect of the Yukawa couplings depends on ln(y r (0)/y°(0)), ln(y&(0)/y£(0)) and ln(y t (0)/y°(0)). 
The contribution of these terms is negative and hence aj,{M z ) is further increased from its MSSM value. 
The reason for this is that the running of Yukawa couplings is similar in both models q], but in the 
extended model they decrease faster due to the fact the gauge couplings in EMSSM at any scale are 
increased from their corresponding MSSM values. This means that, for a fixed low energy input for 
Yukawa couplings and SU(2) and U(l) gauge couplings, the net effect due the presence of the extra 
multiplets is to decrease the Yukawa couplings at unification scale while increasing the value of the 
unified coupling. Hence y(0) < y°(0) for top, bottom and tau quarks and as(M z ) is increased. Note 
that this is true in the case we have the same superpotential in both models which means that the 
coefficients in front of Yukawa couplings in their one loop running equations are the same to those of 
MSSM. 

The Yukawa effects come in eq.(|30|) as log's of ratios y\(0)/y\(Q) with positive sign in front; this 
might no longer be true if the superpotential of the model we consider for the extended model is changed. 
Still, it seems that a domain with large initial Yukawa couplings y\(0) would favour a decreasing effect 
on the strong coupling at electroweak scale, particularly when they are larger than their electroweak 
scale values. This would in turn favour a fast rate of approach towards infrared fixed points for ratios 
of Yukawa to any gauge coupling ||. However, getting a large initial Yukawa coupling, with same low 
energy input values as in MSSM is not easy, due to, as we mentioned it above, the gauge couplings 
competing effect which seems to be difficult to avoid. In cases with unified coupling larger than in 
MSSM, as it happens in our model with and due to the extra heavy states, the gauge couplings increase 
while increasing the scale, causing the Yukawa couplings to decrease as we increase the scale. This could 
eventually be compensated for by increasing the values of the coefficients in front of Yukawa couplings 
in the one-loop running of Yukawa couplings. This is possible if the Yukawa structure is richer than in 
MSSM; such a particular case will be explored in the next subsection. 



7 i.e.they have the same anomalous dimension coefficients 
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The unification scale is given by (see eq.(|25|)) 
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with A = {r, b, t} in this order. The presence of the factor 
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in the expression for M g /M° prevents us from simplifying any further this analytical expression. To 
make predictions one needs, just as in the previous case, the initial value (at unification scale) of the 
Yukawa couplings and their value at the electroweak scale. Since the latter correspond to the masses of 
the third family it is reasonable to take them to be the same in the MSSM and in the extended model. 
Hence K depends on yx(0)/y x (0) only. This dependence is very mild however because the powers g x 
are very small and therefore, with a good approximation we have that 



K « 1 



(36) 



We checked this numerically! 15 ! by using the one-loop running for Yukawa couplings with one-loop running 
for the gauge couplings for MSSM case and for the extended model (in the presence of the extra-matter) 
and observed that the ratios y x (0)/y x (0) are of order 10 giving K = 1 — 2, with larger K (up to 1.8) 
for larger unified coupling. The main contribution to increasing the value of K comes from the bottom 
quark as K contains a negative power of yb(0)/y%(0) < 1, giving an increase factor larger than unity, but 
still very small due to small g x , as mentioned (this is further suppressed by the tau' contribution). For 
n = 1,2,3 case the same discussion about the validity of these results as that presented in jjj applies. 
Hence we have a clear prediction for the unification scale, 
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This increase of the unification scale is very close to that found in reference [Q] and, for a g < 10a° ~ 0.4 
is less than ~ 3.5, and depends on the values of n, with largest value for smallest n. For large n the 
results are stable, as we previously mentioned. We conclude that the effects of Yukawa couplings on 
the scale are small and thus there is no need for a further numerical calculation. 
Finally, the bare mass of heavy states is given by 
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i For generic values for Yukawa couplings for high and low tan/3 case see ref. j^J. 
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Kj = fj + 
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(39) 



(40) 
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with j = {1, 2, 3} and A = {r, b, t}. We again take Yukawa couplings at M z scale be equal in both MSSM 
and our extended model. The dependence of this fi g scale on the Yukawa couplings y\(0)/y1(0) is again 
weak for large n, and therefore the results stay close to those of reference [Q]. We would like to note that 
our results for unification scale, the strong coupling and decoupling scale \x g were all computed in terms 
of Yukawa couplings at the unification and electroweak scale. The values of these Yukawa couplings can 
be determined numerically from their one-loop running only, and we do not need any numerical work 
for the gauge couplings running, thus simplifying the calculation for M g} a^{M z ) and \x g . 



3.3 Family symmetric Yukawa couplings. 

The predictions we make for the strong coupling and the unification scale depend on the type of 
superpotential we assume; this is so because the running of Yukawa couplings (and therefore their 
effect on the running of the gauge couplings) depends on the type of superpotential. We considered in 
the previous section a MSSM type of superpotential and its predictions. We consider now a different 
superpotential, to show how these predictions change. The pattern of the results we present could prove 
helpful in designing models with better phenomenological predictions. 

We thus turn now to a discussion of the implications of a model designed to give an acceptable 
structure for all quark and lepton masses including the light generations 0. In this case the new 
superpotential has the following form: 
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where the structure in the light quark mass matrix is driven via mixing in the Higgs sector so that the 
two light Higgs doublets of the MSSM are mixtures of the Higgs doublets and of fl^ • This form 
is able to reproduce the values of the masses of the third generation for the case the Yukawa couplings 
are given by their infra-red fixed point values in terms of the gauge couplings ||. We use this as our 
starting point and take the fixed point values for the hfj. This means they are all equal, independent 
of Similarly we take hfj and h[j to be i,j independent and at their fixed points. With this we find 



lnZ^(t) 
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where B V {J-) is presented in the Appendix. From eq.(f42[) we obtain the following relations, valid at any 
scale M. larger than fj, g , M. > fj, g 
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For the case Mz < M. < \i g the above type of superpotential is no longer valid and a MSSM-like 
superpotential applies. This is so because the extra heavy states we consider, among which are all Higgs 
fields^ of eq. (41) decouple at [i g . Below this scale fi g , only the third generation' Yukawa couplings 
give important contributions to wavefunction renormalisation coefficients ZY , while for the first two 
generations the factors Z^., j = {1,2} evolve only through their gauge contributions, Z^.,j = {1,2}. 
We thus have that for Mz < M < \i g 
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which, after integration below fi g , gives that 
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Note also the similarity with eq. (pl|), as expected, with the difference that the initial condition for any 
Z Y takes place at [i g scale (with Z Y (/ig) ^ l!) rather than at unification scale. 

Following the procedure outlined above we obtain, in this case, the following change to the strong 
coupling: 
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(46) 



(47) 
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(M z ). Setting all Z equal to unity in (Eq) gives the previous result 



and 5a 3 1 (M z ) = 1 {M Z ) - a c 3 
jlj] where only gauge effects on a^{M z ) were considered, leading to an increase of its value from the 
MSSM value. We also see that the Yukawa effects give in this case two negative signs in eq.(^6|) in the 
square bracket which make its contribution rather small and negative; For comparison see the result 
of eq.(|22]) for the MSSM superpotential. This contribution, if the signs were positive was expected to 
be the largest as Z factors have a steeper running between M g (where they are equal to unity) and 



This implies that the number of extra heavy states is larger than 5 
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\x g than between \i g and M z due to larger coefficients in fl42|) than in (|44|). The second-last and last 
square brackets in eq.(|4"6|) give a small correction to alpha strong and in fact they largely cancel each 
other. We can further express 5a^ 1 (M z ) in function of Yukawa couplings evaluated at M g , at fi g and 
at electroweak scale. 
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(48) 



(49) 



The terms in fl48| ) which contain a g /a° have an overall increasing effect on the strong coupling driving 
it towards values larger than in MSSM as it can be seen by simply plotting their sum with respect to 
olq. We must therefore analyse the effect Yukawa couplings in eq.(EM) have. Their effect on 03 (M 2 ) 



Ig. 

is not always opposite to that of the gauge couplings. In the expression above we see three types of 
contributions. The first one is due to Yukawa effects between unification scale and (effective) decoupling 
scale jj,g. This is due to the family symmetric Yukawa couplings we chose above this scale. We can 
assume that the values of y\(fi) are larger than y\(n g ); this would be expected given that we considered 
the family symmetric couplings, which is possible if they flow to infrared fixed point values, thing 
favoured by a large initial Yukawap^. The overall effect would therefore be to increase the strong 
coupling, due to the negative signs of the last two terms for tau and top couplings in the first curly 

where the log's 



bracket eq.(48) which dominate the bottom contribution. This is unlike equation 
of Yukawa ratios had positive sign in front, (but were themselves negative). Here they are positive, but 
their sign in front is not positive for all bottom, tau, top couplings. 

The second contribution comes from terms like y\(fj, g ) / y\(M z ) and the third contribution comes 
from MSSM terms, y^(0) / y^(M z ) . These terms largely cancel because the running of Yukawa couplings 
below \x g scale is the same in MSSM and EMSSM and we have same electroweak values for y\(M z ), 
y^{M z ) (these are inputs of our analysis); We also know that Yukawa effects in MSSM, coming in as 
ratios y^(0)/y^(M z ) are small and the scale \i g being in general heavyp 1 ], we conclude that the last two 
curly brackets of eq.(|4^) do not solve the problem of how to reduce the small discrepancy between the 
predicted value for the strong coupling and its experimental value. 

To conclude, it is not easy to reduce the strong coupling and increase the unification scale signif- 
icantly, and it seems that the observation we made after eq.(^) is true in general for models which 
consider only complete additional representations. 



20 Considering y\(0) smaller than yx(fig) is not allowed in this model because a rapid flow to infrared fixed points was 
assumed here and hence generation independent Yukawa couplings are possible for y\(Q)/y\{Q) < 1, with Q < M g 
21 The extra states are not protected by any chiral symmetry and therefore their mass is heavy 
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4 The large n limit. 



Our analysis has so far assumed that perturbation theory works well up to the unification scale, and 
that the presence of a given number of states does not affect its convergence. However, for a large 
number of states and a large unified coupling one should carefully consider the limits of the two-loop 
perturbative expansion we applied in this work as well as in ref. Q. It is the purpose of this section to 
explore the phenomenological implications for the case when these limits are reached, in the absence of 
Yukawa effects. 

In the case of large n and large a g the perturbative expansion breaks down. An estimate of where 
breakdown occurs can be obtained by comparing two-loop beta function terms with three loop termsPI 
or, equivalently, one-loop terms in the expansion of the anomalous dimensions of the fields with two-loop 
terms. Although the higher order terms have an additional power of coupling they are also proportional 
to n and in the large n limit this will compensate for the additional coupling. This happens for 
na ~ n 2 a 2 /(47r), or na ~ 0(4ir). In the case of large n the perturbation series in a as well as the 



perturbation series for anomalous dimensions of chiral fields can be resummed [15] to leading order 
0(1/ n). This calculation was done by Jones in [|D| and the reader is referred to this work for full 
details. We will use their results for the resummed anomalous dimensions of the light fields of the 
spectrum for the case the Yukawa couplings are ignored. 

The anomalous dimension for any matter field is given by 
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Here F stands for a matter field in the representation Rp of the gauge group SU (fc)R- The function G 
is defined by 

1 T(3 - 26) sin(7re) 
G ( e ) = ow^ Za2 —— ( 51 ) 
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and has a pole for e = 3/2 which sets the radius of convergence of the resummed series. For G = 1 we 
recover the results of eq. . The running of the gauge couplings is now given by 
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where bij and bj are just the two-loop and one-loop beta functions for the MSSM and where Ai = 0, 
A2 = 2, A3 = 3. /Q are the resummed corrections induced in the large n limit ignoring Yukawa effects 
and are given by 
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The effect of /Cj's on the unification scale and the strong coupling may now be readily computed to give 
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22 For this, expand the denominator of eq. (^) to get three loop terms 

23 For the U(l) factor Ck(Rp) is |F 2 where Y is the usual weak hypercharge jli 
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Table 1: The values of a g , strong coupling at a^(Mz), and M g /M° as a function of (large) n. We 
always have that (na g ) / '(47r) < 3/2. 
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with given in eq.(p3|). To use these equations it is necessary to compute the values of aj(/j, g ). 
We do this by writing the equations equivalent to eq.(52) for a^ l {pL g ) and eliminating the term 
n/(2-7r) In (M g //i g ) between these equations for a^. This gives 
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Note that the quantities JCi depend on ctj(fi g ) as seen from eq. (|53|) . For given a g and large n we solve 



(56) to get aj(fXg) which are then used to compute the ratio M g /M° and a^(M z ) of eqs. (p4|) , (pq). 
The results for values of n and a g which avoid the pole in G are presented in Table 1. 

We see that the effects of a large number of states is very small and the results for M g /M° and 
a^,{M z ) are in general insensitive to the variation of (large) n and stay close to the MSSM values. 



5 The case of strong coupling. 

Our analysis to date has assumed that the couplings (even at large n) remain in the perturbative domain 
between electroweak scale and unification scale. However, non-perturbative physics could be important 
below unification scale and in this section we will consider this possibility. In this case one apparently 
loses all the predictive power because the full beta functions are known only perturbatively and the 
perturbation series does not work beyond the scale of non-perturbative physics. In M it was argued 
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that this is not the case because the ratios of the gauge couplings are driven towards stable infra-red 
fixed points. Since the couplings are initially large, the rate of flow to the infrared fixed point should 
be rapid. As a result, the low-energy values of these ratios are insensitive to their initial values and to 
the non-perturbative effects (provided the couplings initially lie in the domain of attraction of the fixed 
points), giving a reliable prediction for 03 (M z ) even in this case. Using the fixed points (FP) of the 
one-loop RGE to determine the boundary conditions of the gauge couplings at the decoupling scale of 
the additional massive states, it was found that a low value of a^(M z ) is predicted, very close to the 
mean experimental value and in much better agreement than the perturbative MSSM result. 

Here we wish to consider corrections to the above "fixed point" result in order to estimate the 
precision of the result and to establish whether the apparent improvement over the perturbative case 
is significant. To do this we shall continue to use the perturbative solution, but only in the domain 
where it is applicable. As we shall see in this region one may see the RG flow to the fixed points but 
there are calculable corrections. The first point to make is that even using the improved perturbation 
sum discussed in the last section one requires^ na < 0(6-71"). We can only calculate those effects below 
the non-perturbative domain at the scale M Q where the gauge couplings enter the perturbative domain 
and as we have just remarked this starts at quite small values of a. Thus our error estimates will 
necessarily be somewhat rough. There are two corrections to the "fixed point" calculation, deriving 
from the relations among the gauge couplings at the decoupling scale^j fj, g , and from the uncertainty in 
their exact values as they leave the non-perturbative region to enter the perturbative domain. Here we 
analyse them briefly. 

The first correction arises because, even if one sticks to the one-loop beta functions, due to the finite 
energy range involved the ratios of couplings are not driven quite to the fixed point. Therefore the 
boundary condition we used for the evolution of the couplings below fx g is not the "fixed-point" ratio 

Rji(Vg) = \ 9 \ = Rji ( 58 ) 
QiyUg ) 

with = b'j/b'j. Instead, the boundary condition for the running of the gauge couplings below the 
decoupling scalqpl is a "quasi-fixed-point" (QFP) relation for the ratio of the gauge couplings. This 
relation exactly takes into account a one-loop running for the gauge couplings above fj, g scale (just as 
in FP case) and the finite range of energy. This QFP relation is 



Rji(Pg) = 



= (59) 
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This can be easily deduced by integrating the one-loop differential eqs for the gauge couplings above n g . 
Using this QFP boundary condition for the gauge couplings evolution below [i g scale, one gets a value 
for as(M z ) which is typically smaller than that of FP case. The results depend on the value of Rij(M ), 
which enters (|59"|), but this dependence is relatively weak and the prediction for 03 (M z ) stays below 
MSSM value. This is somewhat expected, as we know that two-loop corrections in MSSM increase 
as(M z ) from its rather good one loop prediction^], taking it away from the current experimental upper 
limit. As eq. (59) ignores the two loop evolution' effects between fj, g and M Q , we get a smaller a^(M z ) 



than in MSSM. 



24 We will however ask that na < 0(4n) instead, because we will not use the resummed perturbation series, and this 
relation marks the limit where it becomes important. 

25 In two- loop running of gauge couplings the effective decoupling scale is fig, the bare mass of the heavy spectrum 
26 Below this scale a MSSM like spectrum applies 

27 In MSSM a one-loop calculation gives as{M z ) w 0.117 while two-loop calculation increases it to w 0.126 
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This brings us to the second correction which consists of two- loop or higher order (gauge and Yukawa) 
effects between \x g and M Q . This affects the relation ([59]) and therefore the low energy predictions for the 
strong coupling. We will consider here only the gauge effects. Their effect is to change the expression of 
Rji(fig) such as to increase the value of the strong coupling at electroweak scale. The increase overcomes 
the QFP decrease in a^(M z ), bringing it closer to MSSM prediction. This conclusion is subject to the 
corrections due to the unknown values of the gauge couplings at M D (or equivalently Z(M Q )), where 
they enter the perturbative domain as we lower the scale, but there is a systematic effect increasing the 
value of a^{M z ) even for sizeable changes in the boundary conditions. 

Our analysis is limited because the range between \i g and M Q is quite small and the calculable 
corrections to the fixed point result correspondingly small. However this does not mean the full non- 
perturbative correction to the fixed point result is small. An estimate of these corrections coming from 
the higher order and quasi fixed point terms discussed above may be obtained by noting that there 
is a focussing effect of the fixed point following in the perturbative flow domain which reduces the 
contribution of such effects by the factor 2\m\ ec l(F^I))- Using this factor with ai(M Q ) limited by 
the perturbative condition na < 0(4-71") together with the assumption of 0(1) deviations from the QFP 
boundary condition we arrive at a reasonably conservative estimate for these errors of O(±0.01) for 
n = 12, increasing for larger values of n. 

6 Summary and Conclusions 

For the case of perturbative evolution of the gauge and Yukawa couplings up to the unification scale, 
we derived an analytical method to determine the two-loop unification predictions for the value of 
unification scale, the intermediate mass scale and the value of the strong coupling at M z in models 
with the MSSM spectrum augmented by additional massive representations filling out complete SU(5) 
representations. The effects of the two-loop terms involving Yukawa couplings are in general relatively 
small and model-dependent. For the models we considered, a^{M z ) cannot be lowered below the MSSM 
value, even with Yukawa effects present, keeping az{M z ) above the mean experimental value. However 
the sign of the effect is not universal, so it is possible a richer structure in Yukawa sector could reduce 
the strong coupling. We also showed that the unification scale is not changed significantly by the top, 
bottom and tau Yukawa coupling effects. Using the large n resummed perturbation series we showed 
that the results are stable in the limit of inclusion of a large number of complete SU(5) representations 
and stay close to MSSM predictions. Finally we considered the case that unification occurs at strong 
coupling with perturbative analysis breaking down below the unification scale. Because of the fixed 
point structure in the RGEs (at one- loop level) one may still make a prediction for a^{M z ). However 
threshold and two loop effects above the decoupling scale for the heavy states can be sizeable. As a 
result the low estimate for as (M z ) obtained using fixed point boundary conditions at the decoupling 
scale is increased for large n towards the MSSM value. Overall we estimate an irreducible error in the 
strong coupling determination of as(M z ) of 0(0.01) coming from the residual sensitivity to the initial 
values of the ratios of the couplings as they enter the perturbative domain and from two-loop corrections 
to Rji(fig) above fi g . 

7 Acknowledgments 

D.G. gratefully acknowledges the financial support from the part of University of Oxford and Oriel 
College (University of Oxford). G. A.-C. acknowledges the financial support for his work on this 
project from P.P.A.R.C., the Foundation Blanceflor Boncompagni-Ludovisi and the Swiss National 
Science Foundation. 



17 



8 Appendix 

We have the following expressions for T(R a ) for the fundamental representation 



5 ab T(R) = Tr(T a T b ) = -5' 



arpb\ 



cab 



(60) 



(61) 



For a given flavour the value of T{R a ) is T(R a ) = 1/2. (T(R a ) = 1 for conjugate pairs of fields). 
For the adjoint representation we have 

(T a ) bc = ~if abc 
and therefore, for SU(N) group 

5 ab T(G) = f acd f bcd = N5 ab (62) 
where the structure constants f abc are given by 



abcrpc 



(63) 



The values of A used in text (20) in the context of a MSSM like superpotential and considering only 
the Yukawa couplings for top, bottom and tau is given by (see for example |l6fl) 



A„(F) 



( 

V = T 

V = b 
\ V = t 



h 3 Ql 3 e R 3 U R 3 d R3 Hi H 2 \ 



-1 
-1 
-1 



-2 

-2 



-1 
-2 -3 






"3/ 



(64) 



The running of Yukawa couplings in one-loop, for a MSSM like superpotential can be rewritten, using 
[16|, as follows (we ignore the first two generations' contribution) 



Vb(t) 

V vt(t) J 



I 35 -9 



122 61 



122 



-3 
61 



12 ^2 
61 61 



1 -2 
122 61 



21 
122 



dt 

d\ny b 
dt 



dinyt 
dt 



+ 



305 61 

-23 21 
915 61 



40 \ / «i(t) \ 



61 

160 
183 



V 



136 27 136 
183 



915 61 

For a family symmetric superpotential we used in text (El]) 



a 2 (t) 
V «s(t) J 



+ 0(a 2 ) 



(65) 



we have that 



Bu(F) 











q Lj 


e Rj 


u Rj 


d Rj 


H x 


H 2 


\ 




V 


= r : 


-3 





-6 








-1 









V 


= b : 





-3 








-6 


-3 







V 


V 


= t : 





-3 





-6 








-3 


/ 



(66) 



(67) 



For the above class of superpotential, the Yukawa couplings running in one loop can be written as 
follows 



/ Vr(t) \ 
Vb(t) 

V vt(t) / 



/ 15 -2 1 \ 
/ 146 73 146 \ 



1 16 



-2 20 -5 
219 219 219 



d In y T 
dt 

dt 



\ 



1 -5 13 / \ d\nyt 
438 219 146 / V dt 



^ 73 73 73 \ / ^ 



+ 



7 13 _80_ 

1095 73 219 



OL2{t) 



+ 0(a 2 



(68) 



\ 232 15 232 / \ f + \ I 

\ 3285 73 657 / \ " 3 W / 
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